Cosmic-ray electron injection from the ionization of nuclei 
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We show that the secondary electrons ejected from the ionization of heavy ions can be injected into 
the acceleration process that occurs at supernova remnant shocks. This electron injection mechanism 
works since ions are ionized during the acceleration when they move already with relativistic speed, 
just like ejected electrons do. Using the abundances of heavy nuclei measured in cosmic rays at 
Earth, we estimate the electron/proton ratio at the source to be ~ 10~*, big enough to account for 
the nonthermal synchrotron emission observed in young SNRs. We also show that the ionization 
process can limit the maximum energy that heavy ions can reach. 



Supernova Remnant (SNR) are believed to be the pri- 
mary sources of Cosmic Rays (CR) in the Galaxy. The 
theory of Diffusive Shock Acceleration (DSA) applied to 
SNR blast waves propagating in the interstellar medium 
provides the most comprehensive framework for the ex- 
planation of the CR spectrum measured at Earth. Nev- 
ertheless some aspects of DSA still remain unsolved and 
one of the most fundamental issues concerns how elec- 
trons can be injected into the acceleration process. The 
presence of accelerated electrons is a well established fact, 
deduced from direct observations of young SNRs, where 
both Radio and X-ray emission are interpreted as syn- 
chrotron radiation of highly relativistic electrons [l| . 

The injection problem is related to the shock dynamics 
which is dominated by protons (and maybe also by heav- 
ier ions). The shock layer is expected to be few thermal 
proton gyroradii thick, which means that particles need 
have a few times the mean thermal proton momentum 
in order to cross the shock and undergo DSA. The injec- 
tion condition can be write as p > pinj = ^Pp,th, where 
Pp^th = ■y/2r7ip KbTp^2 is the typical downstream thermal 
proton momentum and ^ contains the complex shock mi- 
crophysics [1| . According to theoretical estimates 0] and 
observational constraints ^ is in the range 2-4. The 
injection condition can be easily fulfilled for suprather- 
mal protons which reside in the highest energy tail of the 
Maxwellian distribution. On the other hand the same 
condition, applied to electrons, is satisfied only in the rel- 
ativistic regime with a minimum Lorentz factor 7inj w 3- 
30 for a typical Tp,2 in the range 10^-10^ K. It is therefore 
hard to imagine how electrons can come from the thermal 
component: if we assume that electrons upstream of the 
shock thermalize downstream their bulk kinetic energy, 
i.e. KBTf,^2 = ^'TieUgjjQj.jj, then the mean thermal elec- 
tron momentum is Pe,th = i'm-e/n^p) Pp,th- Even assum- 
ing some mechanism able to equilibrate quickly electrons 
and protons at the same temperature, the mean electron 
momentum rises only up to Pe,th = ynv/^pPpi'i- 

Most proposed solutions for this injection problem in- 
volve some kind of pre-acceleration mechanism able to 
accelerate electrons from thermal energies up to mildly 
relativistic energies. Some studies predict that electrons 
can be effectively pre-accelerated by electrostatic waves 
generated in the shock layer These mechanisms are 

difficult to study analytically, because their understand- 



ing requires the knowledge of the complex microphysics 
that regulates collisionless shocks. A better way to in- 
vestigate them is through particle-in-cell fi\ and Monte 
Carlo simulations [8], which, unfortunately, are still not 
able to provide firm conclusions on the electron injection 
efficiency 

In this paper we show that the ionization of heavy nu- 
clei during acceleration can inject a number of mildly 
relativistic electrons large enough to account for the syn- 
chrotron radiation observed in young SNRs. In fact, nu- 
clei heavier than hydrogen in the ISM where SNR shocks 
propagate, are never fully ionized simply because the typ- 
ical ISM temperature is not large enough, being of the 
order of 10^ — lO'^ K. This statement is also supported 
by the presence of Balmer-dominated filaments observed 
in several young remnants 0] , showing that even the hy- 
drogen in the ISM is not fully ionized. 

Once the shock encounters a partially ionized atom, 
the atom can start DSA in the same way protons do. We 
note that a correct computation of the injection of heavy 
ions involves the knowledge of the initial charge and the 
aggregate state of atoms and their downstream tempera- 
ture, which are very difficult to predict. We neglect such 
complications and we assume that the injection of heavy 
elements occurs, simply because they are observed in the 
CR spectrum. The key point we want to stress here is 
that, once the acceleration begins, atoms are not stripped 
immediately, because the ionization time turns out to be 
large enough to allow them to reach relativistic energies 
before complete ionization. When atoms move relativis- 
tically, the ionization can occur either via Coulomb col- 
lisions or via photoionization. In both cases, the mean 
kinetic energy of ejected electrons, measured in the ion 
rest frame, is negligible with respect to the electron mass 
energy. Hence ejected electrons move, in the plasma rest 
frame, approximately along the same direction and with 
the same speed of the parent atoms. In this case, the 
momentum of ejected electrons can easily exceed pinj ■ In 
order to prove this statement we start comparing the ac- 
celeration with the ionization time. 

Let us consider a single partially ionized species TV, 
with mean charge Z^ff, atomic charge Z and mass rriN — 
2Zmp. For simplicity we compute the acceleration time 
in the framework of linear shock acceleration theory, i.e. 
neglecting the dynamical role of accelerated particles. 
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and for a plane shock geometry. If a particle with mo- 
mentum p diffuses with a diffusion coefficient D{p\ the 
well known expression for the acceleration time is [10[ : 
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where u is the plasma speed in the shock rest frame, and 
the subscript I (2) refers to the upstream (downstream) 
quantities (note that Ughock = tti)- The downstream and 
upstream plasma speeds are related through the compres- 
sion factor, U2 = ui/r. We limit our considerations to 
strong shocks, which have compression factor r = 4, and 
we assume Bohm diffusion coefficient, i.e. Db = rL^c/S, 
where (3c is the particle speed and = pc/{ZcseB) is 
the Larmor radius. The turbulent magnetic field respon- 
sible for particle diffusion is assumed to be compressed 
downstream according to B2 = r Bi. Even if this relation 
applies only for the magnetic component parallel to the 
shock plane, choosing a different compression rule does 
not affect our main results. Applying previous assump- 
tions, Eq. dl]) becomes: 

taccil) = 1.7 (7 - 7"') u^' {Z/Z,s) yr , (2) 

where 7 is the particle Lorentz factor. Here B^q is the 
upstream magnetic field expressed in pG and the shock 
speed is ui = ug, IO*cm/s. In order to compute the energy 
reached by particles when ionization occurs, we compare 
Eq. ^ with the ionization time. As already mentioned, 
ionization can occur either via Coulomb collisions with 
thermal particles or via photoionization by background 
photons. Whether the former process dominates on the 
latter depends on the ratio between the thermal particle 
and the ionizing photon densities. 

We consider first the role of collisions. A full treatment 
of coUisional ionization is hard because of the complex 
cross section involved and it is beyond the scope of this 
work; our purpose can be achieved using the following 
approximation. Let us consider the process in the rest 
frame of the target atom, which is bombarded with point- 
like charged particles (in our case protons or electrons) 
with kinetic energy -Bkin- The classical ionization cross 
section has a maximum when i?kin is twice the ionization 
energy, /, and the value is: 
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where ao is the Bohr radius, /Ryd is the ionization poten- 
tial expressed in Rydberg units and A^e is the number of 
electrons in the considered atomic orbital. For i?kin > 2/ 
the cross section decreases like E^^^ and reaches a min- 
imum when relativistic effects become important, while, 
in the full relativistic regime tXcoU ck log(i?kin) (see e.g. 
[ll|). Hence we can use Eq. ([3]) as a good upper limit for 
the collisional cross section in a wide range of incident 
particle energy. 

Accelerated ions collide mainly with thermal protons 
and electrons, whose densities are assumed to be equal. 



ne = rip. The total ionization time is the average between 
the upstream and downstream contribution, weighted 
for the respective residence time, i.e. Tcoh = (ti -t- 
t2)(ii/Tcona + i2/Tcoii,2)"\ where = AD^/cui [fO]. The 
final result is: 
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where ni is the upstream proton density in cm^"^. Now, 
equating Tcoh with the acceleration time in Eq. ([2]), we 
get the value of the Lorentz factor, 7coii that ions have 
when the collisional ionization occurs. Using typical pa- 
rameters for a young SNR, the result reads: 
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An upstream magnetic field around 25/iG is expected if 
magnetic amplification occurs, but the condition 7coii > 
7inj can be easily fulfilled even for magnetic field as low 
as the mean galactic value, i.e. 5/iG. 

In the atom rest frame ejected electrons can have 
kinetic energy ranging from up to i?kin — I- But, 
due to the long-range nature of Coulomb interaction, 
events with a small momentum transfer are highly fa- 
vored and the majority of ejected electrons have kinetic 
energy E <C nieC^ [Il[. As a consequence when parent 
atoms move relativistically with respect to the plasma 
rest frame, ejected electrons move approximately in the 
same direction and with the same Lorentz factor com- 
puted in Eq. The same is true when the electron is 
ejected by photoionization. 

Photoionization can occur only when atoms collide 
with photons whose energy is larger than the ionization 
potential /. Atoms moving relativistically see a distribu- 
tion of photons peaked in the forward direction of mo- 
tion, with a mean photon energy e' = 76, where e is the 
photon energy in the plasma rest frame. In order to es- 
timate the photoionization time we adopt the simplest 
approximation for the K-aheW cross section [l^, i.e.: 



tTph(e') = 64a-V^"' (I/e') 
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where ctt is the Thompson cross section and a is the 
fine structure constant. The factor is due to the 
nuclear charge dependence of K-type orbitals dimension. 
In order to get the full photoionization time we need to 
integrate over the total photon energy spectrum, i.e.: 
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where duph/de is the photon spectrum as seen in the 
plasma rest frame. The relevant ionizing photons are 
only those with energy e ~ 1/7 (measured in the plasma 
frame) because the photoionization cross section rapidly 
decreases with increasing photon energy. The corre- 
sponding numerical value of photoionization time is: 



rph(7) ^0.01^2 (nph(//7)/c 



yr. 



(8) 
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Comparing Tph with Eq. ([4]) for the last inner orbital 
(which has Inyd — Z^) we see that photoionization gen- 
erally dominates over the coUisional ionization for heavy 
ions, namely when riph ~ 4niZ~^. 

In Fig. [1] we compare the acceleration time with both 
the coUisional and the photoionization time. The two 
panels show the case of ionization for two hydrogen-like 
ions, He+ and C^+, so that ZcS — Z ~1. Each character- 
istic time is shown for two different choices of the parame- 
ters: iacc, plotted with solid lines, is shown for ui = 3000 
km/s, Bi = 3/iG (upper line) and for ui = lO'* km/s, 
Bi = 20^G (lower line); Tcoh (dot-dashed) is shown for 
Til = 0.01 (upper line) and 1 cm~'^ (lower line); finally 
Tph (dashed line) is computed according to Eq. ([7]), using 
the Galactic interstellar radiation field (ISRF) plus the 
cosmic microwave background. We use the ISRF as cal- 
culated in [l^, which includes the photons produced by 
stars and the infrared radiation resulting from the stellar 
light reprocessed by Galactic dust. In Fig. [T] the lower 
dashed line is computed using the ISRF in the Galactic 
center, while the upper dashed line corresponds to a lo- 
cation in the Galactic plane, 12 kpc far away from the 
Galactic center [l^. We neglect the high energy radia- 
tion coming from the remnant itself because the typical 
number density of the X-ray photons is negligible and 
does not exceed ~ 10~^ph/cm'^. Fig. [1] shows that the 
most relevant contribution to photoionization comes from 
the optical photons, which produces the first dip present 
of the dashed curves. The value of 7 where TcoIi and 
Tph intersect tacc , identifies the Lorentz factor of ejected 
electrons. From the upper panel we see that even elec- 
trons from He+ can be ejected with 7 > 10 and can easily 
undergo DSA. 

A remarkable consequence of the ionization process is 
that, under appropriate circumstances, heavy elements 
reach a maximum energy lower than Z x E^°^°^, the 
value predicted by the shock acceleration theory. This 
because they cannot be completely ionized in a time less 
than the Sedov time of a typical SNS. In fact the pho- 
toionization of the last inner shell, due to photons with 
energy e, occurs only when 7 > Z^Ryd/e. The accelera- 
tion time required to reach such a Lorentz factor is iacc = 

2.4(Z/26)2(e/eV)-i(Bi/20AiG)-i(u/5000km/s)-V- If 
the optical photon density is low enough, the photoion- 
ization is dominated by IR photons (e ^ lO^'^eV), and 
tacc can be longer than the Sedov time. 

Now we estimate whether the number of electrons in- 
jected into the accelerator is large enough to produce 
the observed synchrotron emission. In the literature the 
number of accelerated electrons is usually compared with 
that of protons: DSA operates in the same way for both 
kind of particles, hence a proportionality relation be- 
tween their distribution functions is usually assumed, i.e. 
fe{p) = Kepfp{p) (valid in the energy range where elec- 
tron losses can be neglected). It is worth stressing that 
here we are only interested in the electron/proton ratio 
in young SNRs, and not to the Kf,p measured in the CR 
spectrum at Earth. These two quantities could be dif- 




FIG. 1: Comparison between acceleration (solid lines), pho- 
toionization (dashed) and coUisional ionization time (dot- 
dashed) as functions of the ion's Lorentz factor. The top and 
lower panels show the results for the hydrogen-like ions He^ 
and C^^, respectively. For each time two curves are shown, 
rapresenting two different set of parameters, as explained in 
the text. 



ferent because the latter is the sum of the contribution 
coming from all sources integrated during the source age, 
and also reflects transport to Earth and losses in trans- 
port (expecially radiative losses for electrons). 

The value of Kf,p in the source strongly depends on 
the assumption for the magnetic field strength in the re- 
gion where electrons radiate and, in the context of the 
DSA theory, it can be determined for those SNRs where 
both nonthermal X-ray and TeV radiation are observed. 
Two possible scenarios have been proposed 0, [13, [Tsj . 
In the first one electrons produce both the X-ray and 
the TeV components via the synchrotron emission and 
the inverse Compton effect, respectively; this scenario 
requires a downstream magnetic field around 20/.tG, and 
Kep ^ lO^'^ — 10^"^. The second scenario assumes that 
the number of accelerated protons is large enough to ex- 
plain the TeV emission as due to the decay of neutral 
pions produced in hadronic collisions. In this case the 
DSA requires a magnetic field strength of few hundreds 
^G and K^p ^ 10^^ — 10~^. Such a large magnetic field 
is consistently predicted by the theory as a result of the 
magnetic amplification mechanisms which operate when 
a strong CR current is present. In the following we show 
that injection via ionization can account for the second 
scenario, i.e. the one with efficient CRs production. 
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In order to get the electron spectrum, fe{p), we need 
to solve first the transport equations for all partially ion- 
ized species which take part in the acceleration process 
and release electrons. Then we can use the ions distribu- 
tion functions as source terms for the electron transport 
equation. For the sake of simplicity let us consider the 
simplest case where the acceleration involves only one 
hydrogen-like species which inject one electron per atom, 
iV+ A^++ + e^, as could be the case for He"*". The ac- 
celeration of all the three components can be described 
using the well known transport equation [l^l but addind 
a "decay term" to take into account the ionization pro- 
cess, i.e.: 

dfi ^, .d^fi Idu dfi ^ 

where the index i — e,N^ ^N^^ identifies the species. 
Qi is the source term while Si is the "decay term" due 
to the ionization. We assume that the injection of ions 
occurs only at the shock position and at a fixed mo- 
mentum pinj, hence Qfq+{x,p) — K 6{p~pin-^ 6{x), where 
the normalization constant K is determined by the total 
number of ions injected per time unit. The decay term 
is Sm+ = fN+{x,p)/Tio-a{p), where the total ionization 
time is Tion = {T~Ja + '''ph^)^^- For electrons and 
the decay terms vanish while the injection terms can be 
approximated as follows: 

Q,{x,p) = r dV ^^'^^'f^ 5'^''\p~i^p') ,^^e,N++. 

(10) 

Because both Af++ and e move approximately with the 
same Lorentz factor of N'^, we can set £,i — 1 for Af++ 
and S,i = Trie/m^ for electrons. In the case of linear shock 
acceleration theory, Eq. ([9]) can be solved using standard 
techniques and we will show the detailed procedure 
in a future paper. We define po as the momentum value 
where the ionization time for is comparable to its 
diffusion time yjA^Dju^. It is easy to show that /7V++ {vn) 
and fe{Pe) both become a power law oc for pn > Po 
and pe > Pa x respectively. The index s is only a 



function of the compression factor, s = 3r/{r — 1), and 
s —^ 4: when the shock is strong. In the limit p 3> Po 
the ratio between electrons and ions has the following 
expression: 

- r /e(p) ^ f "^e m\ 
KeN = lim -— = — — - . (11) 

P»Po Jn++{P) 2Z - I V 771 AT / 

Here the factor Z/{2Z — 1) is due to the different diffusion 
coefficient for electrons and ions, while the ratio me/mN 
is due to the different momentum they have when the 
ionization occurs. In order to give an approximated esti- 
mate for Kep we need to multiply Eq. (fTT|) by the total 
number of ejected electrons, i.e. (Z — Zcs), and sum 
over all atomic species present in the accelerator, i.e. 
^ep - J2n ^Np (Zn - Zn^os) KeN, where Knp are the 
abundances of ions measured at the source in the range 
of energy where the ionization occurs. Even if the values 
of are widely unknown, we can estimate it using 

the abundances measured at Earth and adding a correc- 
tion factor to compensate for propagation effects, namely 
the fact that particles with different Z diffuse in a dif- 
ferent way. The diffusion time in the Galaxy is usually 
assumed to be Tdig oc {p/Z)~^, with S « 0.3 — 0.6 (see [1^ 
for a review on recent CR experiments). If i^Afp.o is the 
ion/proton ratio measured at Earth, than the same quan- 
tity measured at the source is K^p = K^p^qZ^^ . Hence 
the final expression for the electron/proton ratio at the 
source is: Kep ~ J2n ^Npfi {Zn - Zn, cs) KeN- 

Using the abundances of nuclei measured at 1 TeV [l3| 
and assuming that Z/2 is the number of electrons effec- 
tively injected by each species, we have Kep ~ 10^*. Re- 
markably this number is exactly the order of magnitude 
required to explain the X-ray emission in the context of 
efficient CR acceleration. We note that the present esti- 
mate is based on linear shock acceleration theory, while 
a correct treatment requires the inclusion of non linear 
effects. 
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